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We present a theoretical investigation of the optical properties of polaronic excitons in stacked 
self-assembled quantum dots, which is based on the non-adiabatic approach. A parallelepiped- 
shaped quantum dot is considered as a model for a self-assembled quantum dot in a stack. The 
exciton-phonon interaction is taken into account for all phonon modes specific for these quantum 
dots (bulk-like, half-space and interface phonons). We show that the coupling between stacked 
quantum dots can lead to a strong enhancement of the optical absorption in the spectral ranges 
characteristic for phonon satellites. 

PACS numbers: 78.67.Hc,73.21.La,73.21.-b 

Non-adiabaticity is an inherent property of exciton-phonon systems in various quantum-dot structures. Non- 
adiabaticity drastically enhances the efficiency of the exciton-phonon interaction. The effects of non-adiabaticity are 
important to interpret the surprisingly high intensities of the phonon 'sidebands' observed in the optical absorption, 
the photoluminescence and the Raman spectra of quantum dots, in particular, an enhancement of these intensities with 
decreasing the quantum-dot size (see, e. g., Refs.[l|0)- Deviations of intensities of the phonon-peak sidebands, observed 
in some experimental optical spectra, from the Franck-Condon progression, which is prescribed by the commonly 
used adiabatic approximation, find a natural explanation within o ur no n-adiabatic approach^*^*^. Recently, stacked 
quantum dots have received increasing attention (see, e. g., Refs. l7l8l l9i tlQlllil2ll3^1 due to the possibility to finely 
control their energy spectra. This makes stacked quantum dots very promising for future nanodevicesi^ii^. In the 
present work, the non-adiabatic approach is applied to stacked InAs/GaAs quantum dots, which reveal a richer 
structure of phonon and exciton spectra in comparison with those for a single quantum dot. 

In order to model coupled self-assembled InAs/GaAs quantum dots, we consider a stack of N parallelepiped-shaped 
quantum dots with sizes Ix, ly, Iz and interdot distance d along the z-axis. Within the present approach, the lateral 
sizes {Ix^ ly) of each quantum dot in a stack are supposed to be much larger than its size Iz along the growth axis. 
The stack is a system of (2A^ + 1) layers (n = 1, . . . 2N + 1) with parameters 

forn = 2,4, ...,2iV; 
forn = 3,5,...,2iV-l; 
for n = l,2iV+ 1. 

The bulk- like optical-phonon frequencies in InAs and GaAs layers of the stacked InAs/GaAs quantum dots coincide 
with the LO-phonon frequencies in InAs and GaAs, respectively. The interface frequencies belong to the stacked 
quantum dots as a whole and satisfy the dispersion equation 

det||afc„(co)|| =0 (fc,n= l,...2iV), (1) 

where akn {'-^) is the dynamic matrix of the interface vibrations with the matrix elements 

ann {(^) = £n (t^) COth + £„+! (uj) coth g|| Z„+i , 
a„,„_i (tj) = a„_i,„ (w) = -£„ (w) /sinhg||?„; 

all other matrix elements are equal to zero. 

In Fig. ^ typical interface-phonon spectra are represented for stacked InAs/GaAs quantum dots formed by two 
InAs parallelepipeds. The frequencies are plotted as a function of the in-plane wave number q||, which takes discrete 
values due to the quantization of the phonons in the xy-plane. In a stack of N quantum dots, each interface-phonon 
frequency of a single quantum dot splits into N branches. The splitting of the interface-phonon frequencies is due to 
the electrostatic interaction between the optical polar vibrations of the different quantum dots. 

These features of the optical-phonon spectrum of stacked quantum dots are manifested in their optical properties. 
We calculate the optical absorption spectrum of polaronic excitons in stacked quantum dots starting from the Kubo 
formula. Within the non-adiabatic approach^ the following expression results for the linear coefficient of the optical 
absorption by the exciton-phonon system in a quantum-dot structure: 
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where f2 is the frequency of the incident Ught, and flp are, respectively, the electric dipole matrix element and the 
Franck-Condon frequency of a transition between the exciton vacuum state and the one-exciton state Exciton 
states in stacked InAs/GaAs quantum dots are determined using an exact diagonalization of the exciton Hamiltonian 
with simple parabolic valence and conduction bands within a finite-dimensional basis of the electron-hole states. The 
evolution operator averaged over the phonon ensemble, U{t), is 
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In Eq. 10}, T is the time ordering 
specific for the quantum-dot structure 
7a (i) are the exciton-phonon interaction 
yx = cxp (-hux/ksT). 

Within the adiabatic approximation, which has been widely used to calculate the optical spectra of quantum dots, 
non-diagonal matrix elements of the exciton-phonon interaction are neglected when calculating a{il) as given by 
Eq. jnj with Eq. Q). In the adiabatic approacbi^*i^ one supposes that (i) both the initial and the final states of a 
quantum transition are non-degenerate, (ii) th e ene rgy differences between the exciton states are much larger than the 
phonon energies. It has been shown in Refs. T'^Wd that these conditions are often violated for optical transitions in 
small quantum dots, which have sizes less than the bulk exciton radius. In other words, the exciton-phonon system in 
a quantum dot can be essentially non- adiabatic. The polaron interaction for an exciton in a degenerate state results 
in internal non-adiabaticity ("the proper Jahn-Teller effect"), while the existence of exciton levels separated by an 
energy cornparable with the LO-phonon energy leads to external non-adiabaticity ("the pseudo Jahn-Teller effect"). 

In Ref. 13, a method was proposed to calculate the absorption spectrum given by Eqs. l|2Jl and 10} taking into 
account the effect of non-adiabaticity on the probabilities of phonon-assisted optical transitions. The key step is 
the calculation of the matrix elements of the evolution operator (/3 |C/(i)| /?'). In order to describe the effect of non- 
adiabaticity both on the intensities and on the positions of the absorption peaks, a diagrammatic approach can be 
used. When calculating these matrix elements we take into account that in a quantum dot, due to the absence of 
momentum conservation, the product (/3i|7aI/32)(/32|7aI/53) can be non-zero for /?i ^ Ps, as distinct from the bulk 
case. Consequently, the evolution operator is in general non-diagonal in the basis of one-exciton wavefunctions |/3). 
For the absorption coefficient we obtain 
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/„(x) is a modified Bessel function of the first kind and Sxp is the Huang- Rhys parameter, which is related to the 
interaction of the exciton in the state P with phonons of the A-th mode: 
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The self-energy terms eIj^-* (Q) and E^^-* (Q) in Eq. Q are obtained by summing diagrams, which describe one- and 
two-phonon non-adiabatic contributions: 

4'^ (^) = E E Pf^p^ " j'^^) M^M^fs (10) 

i=±l A, ft 



and 



^?i^)= E E E ^^;3ft(l^-JlC.Aj 
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The functions Q^pp,{fl) and (3|3^^],(f7) in Eq. (0, which describe contributions of one- and two-phonon processes to 
non-diagonal matrix elements of the evolution operator, take the form 



^ E E (^^ - JUJx) mSftft/S' > (15) 
j = ±l A,ft 
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X (1 - <5/32/3') mSftftft"^S2/33/33/3'-} (16) 



The absorption spectrum is thus expressed through the functions Gfj (Q), which in turn are determined by a closed 
set of equations ©,10, and to (O- 

In Figs. 121 and 13 the calculated optical absorption spectra are shown for a single quantum dot and for a system of 
two stacked quantum dots, respectively. The calculations where performed for low temperatures {{y\ <^ 1}) when the 
absorption-line broadening due to the exciton-LO-phonon interaction is negligible. The broadening shown in Figs. [21 
and 121 is introduced only to enhance visualization. From the comparison of the spectra obtained in the adiabatic 
approximation with those resulting from the non-adiabatic approach, the following effects of non-adiabaticity are 
revealed. First, the polaron shift of the zero-phonon lines with respect to the bare-exciton levels is larger in the 
non-adiabatic approach than in the adiabatic approximation. Second, there is a strong increase of the intensities of 
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the phonon satellites compared to those given by the adiabatic approximation. This increase can be by more than 
two orders of magnitude. Third, in the optical absorption spectra found within the non-adiabatic approach, there 
appear phonon sateUites related to non- active bare exciton states. 

Fourth, the optical-absorption spectra demonstrate the crucial role of non-adiabatic mixing of different exciton and 
phonon states in quantum dots. This results in a rich structure of the absorption spectrum of the exciton-phonon 
systeni^iS*i&. For the stacked quantum dots, this effect is enhanced in the (quasi-) resonant case, when the exciton-level 
splitting, caused by the coupling between quantum dots, is close to a LO phonon energy [see panel (b) in Fig. 
Similar conclusions about the influence of the exciton-phonon interaction on the optical spectra of quantum dots have 
been recently formulated in Ref. ^3 for the "strong coupling regime" for excitons and LO phonons. Such a "strong 
coupling regime" is a particular case of the non-adiabatic mixing related to a (quasi-) resonance, which arises when 
the spacing between exciton levels is close to the LO phonon energy. 

In some cases (see, e. g., Ref. 0) the luminescence spectrum of a quantum dot can be easily derived from its 
absorption spectrum. E. g., under the assumption that the distribution function of the states of an exciton coupled 
to the phonon field, /(fi), depends only on the energy of a state, the luminescence intensity at low temperatures 
{{y\ ^ 1}) can be represented as 



Equation H17|) is applicable, for instance, for thermodynamic equilibrium photoluminescence. In this case the radiative 
lifetime of an exciton is much larger than the time characteristic of radiationless relaxation between one-exciton states. 

We have calculated the spectra of thermodynamic equilibrium luminescence (not shown here) for quantum dots 
with parameters indicated above. Both for single and for coupled quantum dots, the intensities of phonon satellites 
in these spectra are significantly smaller than in the absorption spectrum. This is because in quantum dots under 
consideration the lowest one-exciton energy level (with /3 = 1) is less affected by the phonon- induced non-adiabatic 
mixing of states than higher levels. The luminescence spectrum at low temperatures is dominated by transitions from 
the state with (3 — 1, while the absorption spectrum contains appreciable contributions due to transitions to higher 
exciton-phonon states. 

Due to non-adiabaticity, multiple absorption peaks appear in the spectral ranges characteristic for phonon satellites. 
From the states, which correspond to these peaks, the system can rapidly relax to the lowest emitting state. Therefore, 
in the photoluminescence excitation (PLE) spectra of quantum dots, pronounced peaks can be expected in spectral 
ranges characteristic for phonon satellites. Experimental evidence of the enhanced phonon-assisted absorption due 
to non-adiabaticity has been recently provided by PLE measurements on single self-assembled InAs/GaAsi^ and 
InGaAs/GaAs^^ quantum dots. Our results, which imply significantly more pronounced phonon satellites in PLE 
spectra compared to the luminescence spectra from the lowest one-exciton state, are in line with the experimental 
observationsiSiiS. 
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FIG. 1; Interface-phonon frequencies for two stacked parallelepiped-shaped InAs/GaAs quantum dots. lOq''^^^^ is the frequency 
of LO phonons in InAs at the center of the Brillouin zone. 
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FIG. 2: Absorption spectra, calculated for a single quantum dot with the adiabatic approximation [panel (a)] and with the 
non-adiabatic approach [panel (b)]. Optically active and non-active energy levels of a bare exciton are shown as dotted and 
dashed lines, respectively. Q.\ is the transition frequency for the lowest state of a bare exciton. 
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FIG. 3: Absorption spectra, calculated for a system of two stacked quantum dots with the adiabatic approximation [panel (a)] 
and with the non-adiabatic approach [panel (b)]. Notations are the same as in Fig. 2. 



